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Abstract. We carried out 3D radiative magnetohydrodynamic simulations of the 
convective and magnetic structure in the surface layers (uppermost part of the convec- 
tion zone and photosphere) of main-sequence stars of spectral types F3 to M2. The 
simulation results were analyzed in terms of sizes and properties of the convection cells 
(granules) and magnetic flux concentrations as well as velocity, pressure, density, and 
temperature profiles. 

Our numerical simulations show for the first time a qualitative difference in the mag- 
neto-convection between solar-like stars and M dwarfs. Owing to higher surface grav- 
ity, lower opacity (resulting in higher density at optical depth unity), and more stable 
downflows, small-scale magnetic structures concentrate into pore-like configurations of 
reduced intensity. This implies that in very cool stars magnetic surface structures like 
plage regions and starspots significantly differ from the solar example. Such a differ- 
ence would have major impact on the interpretation of Doppler imaging data and the 
analysis of M dwarf spectra. 



1. Introduction 

Cool main-sequence stars of spectral types F through L have a thick convective enve- 
lope or are fully convective. In many of such stars, magnetic fields of various strengths 
have been detected. In the Sun, the surface magnetic field is observed to be highly 
structured owing to its interaction with the convective flows. This has significant im- 
pact on the magnetic signatures in spectral lines that are used to detect and measure the 
field. In contrast to the Sun, the structure and properties of magnetic fields on other cool 
stars are unknown. In the absence of spatially resolved observations, the effect of the 
magnetic structure on signatures of the magnetic field can be evaluated by numerical 
simulations of the magneto-convective processes. 

2. The MURaM code and simulation setup 

The MHD code MURaM is a "box-in-the-star" code that solves the equations of (non- 
ideal) MHD in three spatial dimensions with constant gravitational acceleration. It 
includes the relevant physical processes such as compressibility, partial ionization, and 
non-grey radiative energy transport. Hyperdiff'usivities and artificial viscosities are in- 
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Figure 1 . Images of the bolometric intensity computed from snapshots of six mod- 
els with stellar parameters according to spectral types F3V, G2V, KOV, K5V, MOV, 
and M2V without magnetic field. Note the different box scale indicated by the yel- 
low bar. The (RMS) intensity contrast (displayed enhanced) is 20.6%, 15.4%, 8.0%, 
7.1%, 3.5%, and 2.1%, respectively. 
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troduced to account for ener gy dissipation on unresolved length scales. For more de- 
tails, see'Vogl er et all (l2005b . 

To fit the surface conditions for different stellar spectral types, gravity and effective 
temperature were adjusted and the opacity bins were recalculated. The size of the com- 
putational box and the spatial resolution were modified in order to cover the relevant 
length scales of the different convection patterns. 

The model grid comprises stars of spectral types F3V (Jeff = 6780 K, logg[cgs] = 
4.301), G2V (Sun;reff - 5770 K, log^[cgs] - 4.438), KO (T^ff = 4950 K, log^[cgs] - 
4.609), K5 (reff = 4370 K, logg[cgs] = 4.699), MO {T,ff = 3910 K, logg[cgs] = 
4.826), and M2 (T^ff = 3690 K, log^[cgs] = 4.826) with solar chemical composition. 
The start models were run with B = until a quasi-stationary state was reached. Then, 
a homogeneous vertical magnetic field with the field strength Bq = 20 G, 100 G, or 
500 G was introduced. 



3. Results 

3.1. General granulation characteristics 

Fig.[T]shows intensity images of snapshots from the six models without magnetic field. 
The typical size of the granules scales roughly with the pressure scale height and thus 
decreases from more than 5 Mm for the F3V model to less than 0.5 Mm for M2V. The 
bolometric intensity contrast decreases from 20% (F3V) to about 2% (M2V). 
As illustrated in Fig.[2l the bright granules are upiiows which carry hotter (brighter) 
gas to the optical surface, whereas the darker lanes between the granules are downfiow 
regions where the radiatively cooled gas sinks back into the convective envelope. The 
typical vertical velocities decrease from about 7 km s~' (F3V) to about 1 km s~' (M2V). 



3.2. Differences in the magneto-convection between G-type and M-type dwarfs 

The modelled magneto-convection shows significant differences between M-dwarfs and 
stars of earlier spectral types, e. g. the Sun. As illustrated by Fig.[3l the initially homo- 
geneous magnetic flux is accumulated into very few structures of high field strength, 
the cause of which are stable downifows. 

While solar magnetic structures appear as bright features, the magnetic structures on 
M-dwarfs tend to be rather dark. In the case of the Sun, magnetic structures create a 
strong depression of the optical surface with hot side walls that can radiatively heat 
the interior of the magnetic structure. This mechanism is much less efficient for the 
magnetic structures in M-dwarf atmospheres: owing to the much higher densities, the 
depressions formed by the magnetic structures are very shallow and their side walls 
have small excess temperatures due to a shallow temperature gradient. Since the mag- 
netic field suppresses convective energy transport, the structures cool down. As Fig. |4] 
shows, magnetic structures in our M2V-star model tend to be darker than the environ- 
ment at all initial fields Bq, the effect being much more pronounced for strong plage 
fields (Bo = 500 G). 

This indicates that plage regions on M-stars might not show bright points but rather 
"pores" and small "star spots" of reduced intensity, which has a crucial impact on the 
interpretation of observational data such as M-dwarf spectra. 
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Figure 2. Vertical velocities v,(rR » 1) at the optical surface computed from a 
snapshot of the M2V-star model (left panel; same time step as in Fig.lTJ. The right 
panel shows the distribution of bolometric intensity /boi(A' = 1) as a function of 
t^v(TR - 1) in a contour plot: contour lines indicate isolines of point density in the 
scatter plot and are drawn for 0.01, 0.03, 0.1, 0.3, 0.5, 0.7, 0.9, 0.97, and 0.99 times 
the maximum point density. 
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Figure 3. Snapshots of the bolometric intensity (left panels) and the vertical com- 
ponent of the magnetic field at the optical surface of the G2V and M2V star models. 
The initial field was Bo = 100 G. The bolometric intensity contrast (displayed en- 
hanced) is 14.9% and 2.30%. 
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Figure 4. Snapshots of the bolometric intensity of the M2V model with different 
initial fields Bq = 20, 100, 500 G (left panels). The intensity contrast (displayed en- 
hanced) is 2.11%, 2.37%, and 5.26% (from top to bottom). The right panels show 
the according distribution of bolometric intensity /boKy" = l)/(^boi(jU - 1)) as a func- 
tion of field strength \B(t x 1)|. Contour lines indicate isolines of the point densities 
in the scatter plot and are drawn for 0.01, 0.03, 0.1, 0.3, 0.5, 0.7, 0.9, 0.97, and 0.99 
times the maximum point density. 



